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Abstract. In this paper, we describe a relation between a microscopic traffic cellular automaton (TCA) 
model (i.e., the stochastic TCA model of Nagel and Schreckenberg) and the macroscopic first-order hydro- 
dynamic model of Lighthill, Whitham, and Richards (LWR). The innovative aspect of our approach, is that 
we explicitly derive the LWR's fundamental diagram directly from the STCA's rule set, by assuming a sta- 
tionarity condition that converts the STCA's rules into a set of linear inequalities. In turn, these constraints 
define the shape of the fundamental diagram, which is then specified to the LWR model. Application of 
our methodology to a simulation case study, allows us to compare the tempo-spatial behavior of both 
models. Our results indicate that, in the presence of noise, the capacity flows in the derived fundamental 
diagram are overestimations of those of the STCA model. Directly specifying the STCA's capacity flows 
to the LWR fundamental diagram, effectively remedies most of the mismatches between both approaches. 
Our methodology sees the STCA complementary to the LWR model and vice versa, so the results can be 
of great assistance when interpreting the traffic dynamics in both models. Especially appealing, is the fact 
that the STCA can visualize the higher-order characteristics of traffic stream dynamics, e.g., the fans of 
rarefaction waves. 

PACS. 45.70.Vn Granular models of complex systems; traffic flow - 47.11.+j Computation methods in 
fluid dynamics - 89.40.-a Transportation 



1 Introduction 

Considering the existing relations between the stochastic 
TCA model of Nagel and Schreckenberg |16U14| and the 
macroscopic first-order model of Lighthill, Whitham, and 
Richards (LWR) jlOU19| . from traffic flow theory, there 
are already numerous links between both modeling ap- 
proaches. An example is the so-called totally asymmetric 
simple exclusion process (TASEP) 4 , which corresponds 
to the LWR model with a noisy and diffusive conserva- 
tion law if a random sequential update is assumed |141 
I15| . Another example is the STCA which can be approx- 



imated by a so-called mean field theory (MFT), and its 
successive refinements, such as the car-oriented mean-field 
theory (COMF), and the recently developed site-oriented 
cluster-theoretic approach 25 ,23 ,24 ,20 21 1 22 . A more 
elaborate discussion on these analytical treatments can be 
found in |llj . In summary, we can say that there already 
exist several methods for bridging both the microscopic 
STCA and the macroscopic LWR model (note that we do 
not consider the class of hybrid models, as we are only 
interested in direct analogies between both microscopic 
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and macroscopic models, and not in pure combinations of 
these model classes). 

In this paper, we reconsider the STCA and LWR mod- 
els, but we take a different approach at studying their re- 
lation: we consider a TCA model as a particle-based dis- 
cretization scheme for macroscopic traffic flow models. It 
is from this latter point of view that our work addresses 
the common structure between both models. Our main 
goal is therefore to provide a means for implicitly incor- 
porating the STCA's stochasticity into the LWR model, 
which is in fact deterministic in nature Note that we 
use the term implicit to denote the fact that the STCA's 
stochasticity is not introduced in the equations by means 
of explicit noise terms. Rather, our methodology implies 
that the stochasticity is introduced through the shape of 
the LWR's fundamental diagram. 

This paper is organized as follows: in section we 
briefly recapitulate both the approaches taken by traf- 
fic cellular automata models (with the STCA in particu- 
lar), and the macroscopic first-order hydrodynamic LWR 
model. In section [31 we then consider a methodology for 
implicitly incorporating the STCA's stochasticity into the 
LWR's triangular fundamental diagram. Continuing, we 
apply this technique to a small case study in section 01 
which points us to some discrepancies between both mod- 
eling approaches. Highlighting some of the resulting arti- 
facts, and investigating the main reason for the difference 
in behavior, we move on to sectional where we present an 
alternate derivation of the fundamental diagrams. Finally, 
the paper concludes with section [HI stating a summary of 
our findings. 

2 Recapitulating the STCA and LWR models 

With respect to the modeling of traffic flows, there are 
largely two model classes possible, i.e., the microscopic 
and macroscopic approach, respectively. In the former class, 
interactions between vehicles in a traffic stream are ex- 
plicitly modeled, giving rise to car-following and lane- 
changing submodels. With respect to the latter class, traf- 
fic streams are mostly treated as inviscid but compress- 
ible fluids. In this section, we briefly recapitulate a special 
class of microscopic models, i.e., traffic cellular automata 
(TCA) models. We then describe a stochastic TCA model, 
called the stochastic traffic cellular automaton (STCA) of 
Nagel and Schreckenberg, after which we conclude with 
an overview of the most prominent features of the macro- 
scopic first-order hydrodynamic model of Lighthill, Whit- 
ham, and Richards (LWR). 

2.1 Traffic cellular automata (TCA) models 

In the field of traffic flow modeling, microscopic traffic 
simulation has always been regarded as a time consum- 
ing, complex process involving detailed models that de- 
scribe the behavior of individual vehicles. Approximately 
a decade ago, however, new microscopic models were being 



developed, based on the cellular automata programming 
paradigm from statistical physics. Let us first describe the 
operation of a single-lane traffic cellular automaton as de- 
picted in Fig. ^ We assume N vehicles are driving on a 
circular lattice containing K cells, i.e., periodic bound- 
ary conditions (each cell can be occupied by at most one 
vehicle at a time). Time and space are discretized, with 
AT = 1 s and AX = 7.5 m, leading to a velocity dis- 
cretization of AV = 27 km/h. Furthermore, the velocity 
Vi of a vehicle i is constrained to an integer in the range 
{0, . . . , fmax}, with Umax typically 5 cells/s (correspond- 
ing to 135 km/h). 
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Fig. 1. Schematic diagram of the operation of a single-lane 
traffic cellular automaton (TCA); here, the time axis is oriented 
downwards, the space axis extends to the right. The TCA's 
configuration is shown for two consecutive time steps t and 
t+ 1, during which two vehicles i and j propagate through the 
lattice. Without loss of generality, we denote the number of 
empty cells in front of vehicle i as its space gap g Si . 

Each vehicle i has a space headway h Si and a time 
headway ht i , defined as follows: 

h Si = g Si + k, (1) 
hi = 9u + Pi- (2) 

In these definitions, g Si and denote the space and 
time gaps respectively; k is the length of a vehicle and pi is 
the occupancy time of the vehicle (i.e., the time it 'spends' 
in one cell). Note that in a traffic cellular automaton the 
space headway of a vehicle is always an integer number, 
representing a multiple of the spatial discretization AX 
in real world measurement units. So in a jam, it is taken 
to be equal to the space the vehicle occupies, i.e., h Si = 1 
cell. 

Local interactions between individual vehicles in a traf- 
fic stream are modeled by means of a rule set. In this 
paper, we assume that all vehicles have the same physical 
characteristics. The system's state is changed through syn- 
chronous position updates of all the vehicles, based on a 
rule set that reflects the car-following behavior. Note that 
most rule sets of TCA models do not use the space head- 
way h Si or the space gap g Si , but are instead based on the 
number of empty cells di in front of a vehicle i. Keeping 
equation Q in mind, we therefore adopt the convention 
that, for a vehicle i its length U — 1 cell; the resulting TCA 
models are called single-cell models. This means that when 
the vehicle is residing in a compact jam, its space head- 
way h Si — 1 cell and its space gap is consequently g Si = 
cells. This abstraction gives us a rigorous justification to 
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formulate the TCA's update rules more intuitively using 
space gaps. 

2.2 The stochastic traffic cellular automaton (STCA) 
model of Nagel and Schreckenberg 

In 1992, Nagel and Schreckenberg proposed a TCA model 
that was able to reproduce several characteristics of real- 
life traffic flows, e.g., the spontaneous emergence of traffic 
jams [TBim] . Their model is called the NaSch TCA, but 
is more commonly known as the stochastic traffic cellu- 
lar automaton (STCA). It explicitly includes a stochastic 
noise term in one of its rules. The STCA then comprises 
the following three rules (note that in Nagel and Schreck- 
enberg's original formulation, they decoupled acceleration 
and braking, resulting in four rules): 

Rl: acceleration and braking 

Vi(t) i- minjv^i - 1) + l,g Si (t - 1), v m ax}, (3) 
R2: randomization 

£(i) < p =► Vi(t) «- max{0, «,•(*) - 1}, (4) 
R3: vehicle movement 

xt(t) «-x<(t-l) +«,•(*). (5) 

The STCA contains a rule for increasing the speed of a 
vehicle and braking to avoid collisions, i.e., rule Rl, equa- 
tion J2J). It furthermore also contains rule R2, equation (J3J, 
which introduces stochasticity in the system. At each dis- 
crete time step t, a random number 6 [0, 1[ is drawn 
from a uniform distribution. This number is then com- 
pared with a stochastic noise parameter p £ [0,1] (called 
the slowdown probability); as a result, there is a probabil- 
ity of p that a vehicle will slow down to Vi(t) — 1 cells/time 
step. According to Nagel and Schreckenberg, the random- 
ization of rule R2 captures natural speed fluctuations due 
to human behavior or varying external conditions. The 
rule introduces overreactions of drivers when braking, pro- 
viding the key to the formation of spontaneously emerging 
jams. Finally, rule R3, equation J5J, allows for the actual 
movement of vehicles in the system. The STCA model is 
called a minimal model, in the sense that all these rules 
are a necessity for mimicking the basic features of real-life 
traffic flows. 

To get an intuitive feeling for the STCA's system dy- 
namics, we have provided two time-space diagrams in Fig-EI 
Both diagrams show the evolution for a global density of 
k = 0.2 vehicles/cell, but with p set to 0.1 for the left 
diagram, and p = 0.5 for the right diagram. As can be 
seen in both diagrams, the randomization in the model 
gives rise to many unstable artificial phantom mini-jams. 
The downstream fronts of these jams smear out, forming 
unstable interfaces |17| . This is a direct result of the fact 
that the intrinsic noise (as embodied by p) in the STCA 
model is too strong: a jam can always form at any density, 
meaning that breakdown can (and will) occur, even in the 



free-flow traffic regime. For low enough densities however, 
these jams can vanish as they are absorbed by vehicles 
with sufficient space headways, or by new jams in the sys- 
tem (Hj. It has been experimentally shown that below the 
critical density, these jams have finite life times with a 
cut-off that is about 5 x 10 time steps and independent 
of the lattice size. When the critical density is crossed, 
these long-lived jams evolve into jams with an infinite life 
time, i.e., they will survive for an infinitely long time jlril 

rum. 



Fig. 2. Typical time-space diagrams of the STCA model (the 
time and space axes are oriented left to right, and bottom to 
top, respectively). The two shown closed-loop lattices each con- 
tain 300 cells, with a visible period of 580 time steps (each vehi- 
cle is represented as a single colored dot). Both diagrams have 
a global density of k = 0.2 vehicles/cell. Left: the evolution of 
the system for p — 0.1. Right: the evolution of the system, but 
now for p — 0.5. The effects of the randomization rule R2 are 
clearly visible in both diagrams, as there occur many unsta- 
ble artificial phantom mini-jams. Furthermore, the speed w of 
the backward propagating kinematic waves decreases with an 
increasing p. 



2.3 The macroscopic first-order traffic flow model of 
Lighthill, Whitham, and Richards (LWR) 

Considering traffic as an inviscid but compressible fluid, 
it follows from this assumption, that densities k, space- 
mean speeds v s , and flows q are defined as continuous 
variables, in each point in time and space, hence leading 
to the names of continuum models, fluid-dynamic models, 
or macroscopic models. 

The first aspect of such a fluid-dynamic description of 
traffic flow, consists of a scalar conservation law; 'scalar' 
because it is a first-order partial differential equation (PDE) 
A typical derivation can be found in [5] and : the deriva- 
tion is based on considering a road segment with a finite 
length on which no vehicles appear or disappear other 
than the ones that enter and exit it. After taking the in- 
finitesimal limit, this will result in a PDE that expresses 
the interplay between continuous densities and flows on a 
local scale: 

dk(t,x) dq(t,x) _ 

3t + dx ' 1 ' 

with the density k and flow q dynamically (i.e., time 
varying) defined over a single spatial dimension. Lighthill 
and Whitham were among the first to develop such a traf- 
fic flow model in 1955 JUj. One year later, Richards in- 
dependently derived the same fluid-dynamic model |19| . 
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albeit in a slightly different form. Because of the nearly 
simultaneous and independent development of the theory, 
the model has become known as the LWR model, after 
the initials of its inventors who receive the credit. In some 
texts, the model is also referred to as the hydrodynamic 
model, or the kinematic wave model (KWM), attributed 
to the fact that the model's solution is based on char- 
acteristics, which are called kinematic waves (e.g., shock 
waves) . 

Crucial to their approach, was a fundamental hypoth- 
esis, essentially stating that flow is a function of den- 
sity, i.e., there exists a q e (k(t, x)) equilibrium relationship, 
more commonly known as the fundamental diagram [Jj]. 
Central to their theory, Lighthill and Whitham assumed 
that the fundamental hypothesis holds at all traffic densi- 
ties, not just for light-density traffic but also for congested 
traffic conditions. 

In order to solve equation @ , we also need the funda- 
mental relation of traffic flow theory, which relates the 
macroscopic traffic flow characteristics density k (vehi- 
cles/kilometre), space-mean speed v s (kilometres/hour), 
and flow q (vehicles/hour) to each other as follows: 



q = k v s . 



(7) 



In general however, there are two restrictions, i.e., the 
relation is only valid for (1) continuous variables, or smooth 
approximations of them, and (2) traffic composed of sub- 
streams (e.g., slow and fast vehicles) which comply to the 
following two assumptions: 

Homogeneous traffic 

There is a homogeneous composition of the traf- 
fic substream (i.e., the same type of vehicles). 

Stationary traffic 

When observing the traffic substream at dif- 
ferent times and locations, it 'looks the same'. 
Putting it a bit more quantitatively, all the vehi- 
cles' trajectories in a time-space diagram should 
be parallel and equidistant 

The latter of the above two conditions, is also referred 
to as traffic operating in a steady state or at equilibrium. 
Employing the fundamental diagram q e (k), relates the two 
dependent variables in equation to each other, thereby 
making it possible to solve the partial differential equa- 
tion. Thus, reconsidering equation JJJJ, taking into account 
the fundamental diagram, the conservation law is now ex- 
pressed as: 



h + q e (k)x = 0, 



(8) 



in which we introduced the standard differential cal- 
culus notation for PDEs. Recognizing the fundamental re- 
lation of traffic flow theory 10, the conservation law JSJ 
can also be cast in a non-linear wave equation, using the 
chain rule for differentiation pj: 



ak 



Analytically solving the previous equation using the 
method of characteristics, results in shock waves that travel 
with speeds equal to: 



dg e (k) 
dk 



(10) 



i.e., the tangent to the q e {k) fundamental diagram. 
As a consequence, solutions, being the characteristics, of 
equation have the following form: 



k(t, x) = k(x — wt), 



(11) 



with the observation that the density is constant along 
such a characteristic. Whenever in the solution of the 
conservation equation, two of its characteristics intersect, 
the density takes on two different values (each one be- 
longing to a single characteristic). As this mathematical 
quirk is physically impossible, the so-called entropy solu- 
tion states that both characteristics terminate and breed 
a shock wave; as such, these shock waves form boundaries 
that discontinuously separate densities, flows, and space- 
mean speeds 5 . The speed of such a shock wave is related 
to the following ratio |T5) : 



w shock 



Aq 

Ak' 



(12) 



0. 



(9) 



with Aq = q u — q^ and Ak = k u — kd the relative differ- 
ence in flows, respectively densities, up- and downstream 
of the shock wave. 

Note that going from a low to a high density regime 
typically results in a shock wave, whereas the reverse tran- 
sition is accompanied by an emanation of a fan of charac- 
teristics (also called expansion, acceleration, or rarefaction 
waves). In shock wave theory, the densities on either side 
of a shock are well defined (i.e., unique solutions exist); 
along the shock wave however, the density jumps discon- 
tinuously from one value to another. 

Because of the well-defined properties of the LWR model, 
it is possible to derive analytical solutions to certain types 
of problems. These solutions can even be drawn graphi- 
cally in a time-space diagram, thereby clearly showing the 
evolution of first-order macroscopic traffic flow character- 
istics (e.g., the speed of backward-travelling jams, ...). 
Besides the previous analytic derivation of a solution to 
the conservation law expressed as a PDE, it is also possi- 
ble to treat the problem numerically. Converting the PDEs 
into finite difference equations (FDEs), and solving them 
numerically stable, can be done by casting the LWR model 
in the context of Godunov FDE methods, allowing for arbi- 
trary q e (k) fundamental diagrams - As such, the LWR 
model is very appealing because it provides an attractive 
means for gaining insight into the primary characteristics 
of the evolution of traffic flows. As a surplus, the model has 
the benefits of the ability to efficiently construct analytical 
solutions, as well as the existence of thoroughly-developed 
and understood numerical schemes. 
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3 Implicitly incorporating the STCA's 
stochasticity 

As mentioned in the introduction, we reconsider the STCA 
and LWR models, taking a different approach at studying 
their relation. Our main goal is to provide a means for 
implicitly incorporating the STCA's stochasticity into the 
LWR model. To this end, we provide a practical method- 
ology for specifying the fundamental diagram to the LWR 
model. Assuming that a stationarity condition holds on 
the STCA's rules, we incorporate the STCA's stochastic- 
ity directly into the LWR's fundamental diagram. 

Relating both the STCA and the LWR models is now 
done using a simple two-step approach, in which we first 
rewrite the STCA's rules into a single rule, leading to a set 
of linear inequalities. These constraints can be considered 
as a v Sc (h s ) fundamental diagram (see e.g., the left part 
m Fig.E|). This latter dia gram can then be converted into 
an equivalent triangular flow versus density q e (k) funda- 
mental diagram. 

3.1 Rewriting the STCA's rule set 

Considering a vehicle's average speed, the STCA's rules 
Rl and R2, equations 10 and (0} respectively, state that 
a vehicle slows down with probability p, and that it does 
not slow down with probability 1 — p. As such, they can be 
rewritten into the following single rule that is expressed 
in continuous speeds and space gaps: 

Vi{t) <— p ■ min{v 4 (i - 1) >T ^-T, g Sl (t - 1) - 1, 
Umax - 1} + (1 -p) ■ mm{vi(t - 1) + 1, 
9s t (t - l),-u m ax}, (13) 

with Vi(t) <— max{0, Vi(t)}. Furthermore, the following 
two algebraic relations always hold: 

a ■ min{&, c} = min{a6, ac}, (14) 
min{a, b} + min{c, d} = min{a + c, a + d, 

b + c,b + d}. (15) 

Applying relation l|14|l to our rule l|13|) . yields the fol- 
lowing result: 



p(g Si (t - 1) - 1) + (1 - p)( Vi (t - 1) + 1), 
p(g Si (t-l)-l) + (l-p)g Si (t-l), 
P(9s, (t - 1) - 1) + (1 - p)vm&x, 
p{v m ax - 1) + (1 - p)(vi(t - 1) + 1), 
p(vma,x - 1) + (1 - p)g Sz [t - 1), 
p(vma,x - 1) + (1 - p)v max }. (17) 

Expanding all the terms between parentheses gives the 
following result: 

Vi(t) <- min{ £y4^T}~+ m{t - 1) + 1 - ^jJ^T) - p, 
pvi{t-l)+g si {t-l)-pg si {t-l), 
pvi(t - 1) + Umax - pwmax, 
pg Si (t-l)-p + Vi(t-l) + l- pvi(t -l)-p, 

pg Si (t-l)-p + Umax - JWmax, 
pvm&x -p + Vi(t-l) + l-pVi(t-l)-p, 
pvma,x -p + g Sl (t -l) - pg s% (t - 1), 
pymsx-p + vmax -gU&xtx}- (18) 

And finally, regrouping for p yields: 

Vi(t) <— min{ Vi{t - 1) + 1 — p, 

p(Vi(t-l)-g ai (t-l))+g ai (t-l), 

p{vi(t - 1) - ti max ) + Umax, 

p(g Si (t - 1) - Vi(t - 1) - 2) + Vi(t - 1) + 1, 

-l)~P, 
p{g Si {t - 1) - w max - 1) + Umax, 
p{vm&x - Vi(t - 1) - 2) + Vi(t - 1) + 1, 
p{vmax - g Si (t - 1) - 1) + g 8i (t - 1), 
vma,x-p}- (19) 

If we now assume traffic is stationary (see e.g., Da- 
ganzo's description of stationary traffic in section l2~3*|l . 
then we can assert that the state of a vehicle at time t 
is the same as its state at time t — 1, i.e., Vi{t) = Vi(t — 1) 
and g Si (t) = g Si (t — 1). As a result, equation Ijl9|l gets 
transformed into the following set of linear inequalities 
that express constraints on the relations between Vi(t), 
g Si {t) 7 p, and wmax: 



Vi(t) <- min{pvi(t- l),p(g Si (t- 1) - l),p(v ma ,x - 1)} 
+min{(l - p){ Vi {t - 1) + 1), (1 -p)g Si (t - 1), 
(l-p)wmax}- (16) 

Using relation (|15|l to this result, allows us to obtain 
a formulation with a single minimum-operator: 

Vi {t) «- min{ pv.it - 1) + (1 -p)(Vi(t - 1) + 1), 
pv i {t-l) + (l-p)g Si (t-l), 
pv t (t - 1) + (1 - p)j)max, 



P{vi{t) -9sM) +9sM > Vi{t) (C2), 

p(vi(t) - w max ) + "max > Vi(t) (C3), 

p(g Si {t)-Vi{t)-2)+^+l>^ (C4), 

9sM-P><t) (C5), 

p(9si(t) - "max - 1) + "max > Vi(t) (C6), 

p(«max-«i(i)-2)+jfc^r+l>jte^T (C7), 

p(vm&x-g s ,(t) - l)+g.M > Vi(t) (C8), 

wmax ~P > Vi(t) (C9). 
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Let us now examine each of these nine constraints CI 
through C9. 

• Constraint CI states that 1 — p > 0, i.e., p < 1. This 
logically follows from the STCA's condition that p G 
[0,1]. ' 

• Constraint C2 states that p{vi(t) — g Si (t)) + g Si (t) > 
Vi(t), i.e., g Si (t)£l^-pj > Vi(t)£L^-pJ. This corresponds 
to Vi(t) < g Si (t), which states that vehicles strive for 
collision-free driving. 

• Constraint C3 states that p{vi{t) — "max) + %iax > 
Vi(t), i.e., UmaxXL— > v i{t)J^-^Pj- This corresponds 
to Vi{t) < «max, which logically follows from the STCA's 
condition that Vi{t) 6 {0, . . . , «max}- 

• Constraint C4 states that p(g Si (t) — Vi(t) — 2) + 1 > 0, 

i.e., Vi(t) <<? 5i (t)-2+i (forp/0). 

• Constraint C5 states that g Si (t)—p > t>.;(i), i.e., Vi(t) < 
g Si (t) —p, which is a more stringent constraint than C2 
and C4. 

• Constraint C6 states that p(g Si (t)— v max — 1)+Vmax > 
Vi(t), i.e., Vi{t) < u max (l -p) +p(g Sl (t) - 1). 

• Constraint C7 states that p(vmax — Vi(t) — 2) + 1 > 0, 
i.e., Vi(t) < Umax - 2 + ± (for p ^ 0). 

• Constraint C8 states that p(wmax—5si(0 — l)+5si (^) — 
u;(t), i.e., < g Si (i)(l -p)+p(vraax ~ !)■ 

• Constraint C9 states that umax-p > Vi(t), i-e., «i(t) < 
^max — P, which is a more stringent constraint than 
C3 and C7. 

Taking the previous considerations into account, we 
can see that constraints CI, C2, and C3 are always sat- 
isfied. The remaining three pairs of similar constraints on 
the relations between Vi(t), g Si {i), p, and wmax, are the 
following: constraints C5 and C9, C4 and C7, and C6 and 
C8. 

In order to gain insight into the more difficult con- 
straints C6 and C8, we first rewrite them as follows: 




■Umax + 1 g ai (t) 



«max 



9s At) 



Fig. 3. A visual representation of the constraints C6 and C8. 
Left: as p — — » 1, C6 changes from a horizontal line at Vi(t) = 
Umax, to a slanted line with a slope of +1, intercepting the 
horizontal and vertical axes at +1 and —1, respectively. Right: 
at the same time, constraint C8 changes from a slanted line 
with a slope of +1, passing through the origin, into a horizontal 
line at Vi(i) = Umax — 1. 



(C6) v l (t) < ^p^ g Si (t) + (1 -p)vma,x - P , 

slope intercept 
(C8) Vj(t) < (1 -p)g a M + p(«max - 1) , 
slope intercept 

where we have separated the terms containing g Si (t) . 
Plotting the speed versus the space gap g Si (t) in 

Fig. |3 allows us to more easily interpret the combined 
effects of these two constraints. On the one hand, if we 
continuously change p = — ► 1, then constraint C6 goes 
from a horizontal line at vi(t) — Vmax, to a slanted line 
with a slope of +1, intercepting the horizontal and vertical 
axes at +1 and —1, respectively. In all cases, the point at 
(i>max + Mmax) remains invariant. On the other hand, 
changing p = — > 1 turns constraint C8 from a slanted 
line with a slope of +1, passing through the origin, into a 
horizontal line at Vi{t) — «max — 1. In all cases, the point 
at (i>max — l,«max — 1) remains invariant. 



3.2 Deriving the fundamental diagram 

The next step of our approach, considers the most deter- 
mining linear inequalities C5, C6, C8, and C9 as bound- 
aries in a v Sc ('g s ) fundamental diagram. As such, we note 
the following observations: 

• Increasing the slowdown probability p, holding 
Umax constant: 

— The average speed wg- in the free-flow regime de- 
creases towards Umax — P- 

— The transition point at the critical space gap g Sc 
remains invariant. 

— The space gap g Sj , corresponding to the jam den- 
sity, increases. 

• Decreasing the maximum speed «max 5 holding 
p constant: 

— The average speed in the free-flow regime de- 
creases towards %ax — P- 

— The transition point at the critical space gap g Sc 
decreases. 

— The space gap g Sj , corresponding to the jam den- 
sity, remains invariant. 

From equation JJl from traffic flow theory, it follows 
that a vehicle's space headway h s is equal to its space gap 
g s (i.e., the distance between the vehicle's frontal bumper 
and the one of its direct frontal leader) , plus the vehicle's 
own length /. As such, the derived v Sc (g s ) fundamental 
diagram can be converted into a v Sc (h s ) fundamental dia- 
gram. Because we originally started from a single-cell TCA 
model (i.e., the STCA model with all vehicles having the 
same unit length) , we can use our convention which states 
that a vehicle's length h > 1 cell cx AX (see our discussion 
at the end of section l2~T1 for more details). 

Let us now consider the relation between the macro- 
scopic traffic flow characteristic density k and the mi- 
croscopic characteristic average space headway h s , i.e., 
h s = fc~l |26ll3*] . This allows us to effectively transform 
the v Sc (h s ) fundamental diagram into a v s<! (k) fundamen- 
tal diagram. As can be seen in the left part of Fig. in- 
creasing the stochasticity leads to the same observations 
that we previously mentioned. Finally, using the funda- 
mental relation of traffic flow theory Q (see section 12.31) , 
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our constraints are transformed into an equivalent trian- 
gular q e (k) fundamental diagram. Applying this technique 
results in the following analytical expressions for the pa- 
rameters of the LWR's fundamental diagram: 



vq = (ttaax —p) —rL 3.(1. 



AT 



1000 



1000 



(20) 
(21) 

(22) 



with I = 1 cell as the length of all vehicles in the single- 
cell STCA model. The capacity flow is calculated using 
the fundamental relation JJJ, resulting in the following 
expression: 

<7cap = fc C rit %• ( 23 ) 
As is visible in the right part of Fig. 0] an increase of 
the stochasticity leads to a lower capacity flow qcap, an 
invariant critical density k c , and a smaller jam density kj. 





Fig. 4. Left: deriving a stationary v 3e (h s ) fundamental di- 
agram from the STCA's constraints CI - C9. The stochastic 
diagram has a higher inverse jam density, but the same inverse 
critical density as its deterministic counterpart (for the same 
Umax)- Right: an equivalent triangular q e (k) fundamental dia- 
gram. 

In conclusion, we note how rewriting the STCA's rule 
set allowed us to obtain a stationary triangular q e (k) fun- 
damental diagram. This fundamental diagram, which im- 
plicitly incorporates the STCA's stochasticity, can then 
be specified as a parameter to the macroscopic first-order 
LWR model of section 12.31 Finally note that the shape 
of the derived fundamental diagram is dictated by the 
inequality constraints CI - C9. As such, it actually repre- 
sents an 'outer envelope', that is to say, all possible funda- 
mental diagrams lie beneath this envelope. This includes 
curved fundamental diagrams, more generally piecewise 
linear fundamental diagrams, et cetera. 



4 Application to an illustrative case study 

After deriving a relation between the STCA and LWR 
models by means of the process explained in the previous 
section we now apply our methodology to a small case 
study. We first describe the setup of the test scenario, after 
which we interpret and discuss our obtained results. 



4.1 Description of the case study 



The case study we consider, is modeled as a single-lane 
road that has a middle section with a reduced maximum 
speed (corresponding to e.g., an elevation, a speed limit, 
. . . ) . This road consists of three consecutive segments A, 
B, and C, as depicted in Fig. 03 whereby vehicles enter 
the road at segment A, travel through segment B, and 
exit it at the end of segment C. For the STCA, we assume 
a temporal and spatial discretization of AT —Is and 
AX = 7.5 m, respectively. The first road segment A then 
consists of 1500 cells (11.25 km), while the second and 
third segments B and C each consist of 750 cells (i.e., 
each approximately 5.6 km long). The maximum speed for 

segments A and C is i> m 'ax = 5 cells/time step, whereas 
1 cell/time step for segment B. The capacity 



it is w max 



flows for all three segments are denoted as g^ap and qc & p 



km 



1 1 km 



17 km 



22 km 



traffic 
» 

demand 



^135 km/hj 



( 2Tkm/h ) 



^135 km/hj 



Section A 



Section B 



Section C 



Fig. 5. The single-lane road of the case study we consider, 
consisting of three consecutive segments A, B, and C. As- 
suming temporal and spatial discretizations of AT = 1 s and 
AX = 7.5 m, respectively, segment A is composed of 1500 cells, 
while segments B and C are each composed of 750 cells. The 
maximum speed for segments A and C is Umax = 5 cells/time 
step, whereas it is Umax = 1 cell/time step for segment B. 



This road is simulated using both the STCA and the 
LWR model, each time for 3000 time steps. As for the 
boundary conditions, we assume an overall inflow of ?cap/^' 
except from time step 200 until time step 600, where we 
have created a short traffic burst of increased demand, 
with an inflow of (<?c a p + <Zcap)/2- Fig. El shows a close up 
of the individual vehicle trajectories for the STCA in a 
time-space diagram, near the border between segments A 
and B. As can be seen from the trajectories, heavy con- 
gestion sets in and flows upstream into segment A, where 
it starts to dissolve at the end of the traffic burst. The re- 
sult is a typical triangular-shaped region that contains a 
queue of slow-moving vehicles (the backward propagating 
waves are clearly distinguished as the pattern of parallel 
black and white stripes). 

Applying our previously discussed methodology, we 
construct a stationary triangular q e (k) fundamental di- 
agram. Its parameters are calculated by means of equa- 
tions (J2DJI - Th° results are listed in Table U with 
the TCA's parameters expressed in cells/time step, ve- 
hicles/cell, and vehicles/time step, respectively, and the 
LWR's parameters expressed in kilometres/hour, vehicles/- 
kilometre, and vehicles/hour, respectively. 
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Fig. 6. A close up of the individual vehicle trajectories for 
the STCA in a time-space diagram, near the border between 
segments A (w^ax = 5 cells/time step) and B («max — 1 
cell/time step), for p — 0.1 everywhere in the system. We can 
see the formation and dissolution of an upstream growing con- 
gested region at the end of segment A, related to the short 
traffic burst. 



"max = 1 (p = 0.1) 


Umax = 5 (p = 0.1) 




TCA 


LWR 




TCA 


LWR 


v s 

k ., 

"'crit 

k- 

^jam 
Qcap 


0.90 
0.50 
0.91 
0.45 


24.30 
66.67 
121.20 
1620.08 


v s 

^crit 
k- 

"■jam 
9cap 


4.90 
0.17 
0.91 
0.83 


132.30 
22.22 
121.20 
2939.71 




«max = 1 (p — 0.5) 


«max = 5 (p = 0.5) 




TCA 


LWR 




TCA 


LWR 


v s 
^crit 

k- 

"•jam 
Qcap 


0.50 
0.50 
0.67 
0.25 


13.50 
66.67 
88.89 
900.05 


v s 
^crit 

k- 

"■jam 
9cap 


4.50 
0.17 
0.67 
0.77 


121.50 
22.22 
88.89 
2699.73 



Table 1. The resulting parameters for the triangular funda- 
mental diagrams, as calculated by means of equations I20H - 
12311 . The TCA's parameters are expressed in cells/time step, 
vehicles/cell, and vehicles/time step, respectively, whereas the 
LWR's parameters are expressed in kilometres/hour, vehi- 
cles/kilometre, and vehicles/hour, respectively. 



4.2 Results and discussion 

The result of numerically solving the LWR model for the 
case of p = 0.1 using the Godunov method [2111], is de- 
picted in the right part of Fig. Note that for the LWR 
model, each cell in the Godunov scheme corresponds to 5 
(i.e., Wmax) consecutive cells of the STCA model. Com- 
paring the tempo-spatial behavior of the LWR model to 
that of the microscopic system dynamics of the STCA 
model (i.e., the left part of Fig.EJ), we find a good qualita- 
tive agreement between the two approaches. With respect 
to the first-order traffic flow characteristics, we note that 
the buildup and dissolution of congestion queues are fairly 
analogous for both techniques. 

In Fig. |S1 we show the results when repeating the same 
experiment, but this time with the stochastic noise p set 
to 0.5 for all three segments. As revealed by the shape of 
the dark triangular region in the LWR model (right part), 
the buildup and dissolution of congestion queues seems to 




Fig. 7. Time-space diagrams showing the propagation of den- 
sities during 3000 time steps for the road in the case study. Left: 
the microscopic system dynamics of the STCA model. Right: 
the results for the LWR model. In both cases, p = 0.1, with 
darker regions corresponding to more congested traffic con- 
ditions. There is a qualitatively good agreement between the 
two approaches on the level of first-order traffic flow charac- 
teristics: the buildup and dissolution of congestion queues are 
fairly analogous for both techniques. 

be exaggerated, especially in the upstream flowing queue 
of segment A. 




Fig. 8. Time-space diagrams showing the propagation of den- 
sities during 3000 time steps for the road in the case study. Left: 
the microscopic system dynamics of the STCA model. Right: 
the results for the LWR model. In both cases, p — 0.5, with 
darker regions corresponding to more congested traffic condi- 
tions. As revealed by the shape of the dark triangular region 
in the LWR model (right part), the buildup and dissolution of 
congestion queues is exaggerated, especially in the upstream 
flowing queue of segment A. 

It is interesting to note that the STCA model reveals a 
higher- order effect that is not visible in the LWR model: 
there exists a fan of forward propagating density waves 
in segment B (see the left parts of Fig. and Fig. [SJ. As 
such, in its tempo-spatial diagram, the STCA seems to 
be able to visualize the characteristics that constitute the 
solution of the LWR model. 

In order to more rigourously quantify the discrepancies 
between the time-space diagrams of both STCA and LWR 
models, we provide their absolute differences in Fig. 
The left part shows the differences for p = 0.1, whereas 
the right part shows the differences for p = 0.5. The most 
important features to look at, are the dark colored regions 
which indicate larger differences between both modeling 
approaches. We can clearly see that there is a problem 
with respect to a quantitative agreement between both 
STCA and LWR models. It appears as though the LWR 
model overestimates the STCA's capacity flows. As a re- 
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suit, it dissolves its jams more quickly in segment B, and 
it predicts a more severe onset of congestion in segment 
A (i.e., the triangular-shaped region containing the spill- 
back queue is more pronounced in the LWR's case). The 
sharply pronounced darker regions in the tempo-spatial 
left part of segment B, are due to the fact that the LWR 
model does not visualize the characteristics of its solution, 
in contrast to the STCA model which is able to give a clear 
indication of them. 




noise terms in the LWR equations (e.g., the conservation 
equation), or as a noise term in the q e {k) relation (i.e., 
the fundamental diagram). We refrain from changing the 
LWR's conservation equation, because this amounts to in- 
troducing some form of numerical diffusion, similar to the 
viscosity terms that are originally encountered in the con- 
servation equation's © right-hand side (which we set to 
zero in order to obtain the inviscid LWR model). 

In Fig. If 01 we show the results of supplying uniformly 
distributed additive noise of O.f (left part) and 0.5 (right 
part). As can be seen, the introduction of noise in the 
fundamental diagram, leads to a 'spreading' of the solu- 
tion. For small noise levels, some of the characteristics 
are revealed; for larger noise levels, the characteristics are 
clearly pronounced, including long jam dissolution times. 



Fig. 9. Time-space diagrams showing the differences in densi- 
ties for the STCA and LWR models, during 3000 time steps for 
the road in the case study. Darker regions indicate large differ- 
ences between both modeling approaches. Left: the differences 
for p = O.f are less pronounced, showing only a dark edge at 
the bottom triangular-shaped region in segment A. Right: the 
differences for p = 0.5, showing significant discrepancies in the 
bottom of the triangular-shaped region in segment A. 

One of the main reasons for this discrepancy between 
both modeling approaches, lies in the derivation of a trian- 
gular q e (k) fundamental diagram for the LWR model, as 
was explained in section^ Because we assumed a station- 
arity condition on the STCA's rule set, the resulting con- 
straints implied an invariant critical density, and always 
overestimated the STCA's capacity flows. In our opinion, 
the different behavior of both models, mainly stems from 
this artifact. As a result, the discrepancies will become 
more articulated when increasing the stochastic noise p. 

5 Alternate derivation of the fundamental 
diagram 

Considering the results of the previous approach, i.e., de- 
riving the LWR's fundamental diagram based on the STCA's 
rule set, and the problems related to it, the next step is 
to specify the fundamental diagram directly, based on the 
empirically observed behavior of the STCA model. In the 
following two sections, we first discuss the effects of explic- 
itly adding noise to the LWR's fundamental diagram, after 
which we discuss our obtained results when specifying the 
fundamental diagram directly. 

5.1 The effect of adding noise to the LWR's 
fundamental diagram 

Adding noise to the LWR model can mainly be accom- 
plished via two ways: either by explicitly incorporating 




Fig. 10. Time-space diagrams showing the propagation of 
densities during 3000 time steps for the road in the case study. 
Depicted are the results for the LWR model, with noise levels 
of 0.1 (left) and 0.5 (right). Higher noise levels clearly reveal 
the typical characteristics of the solution, and introduce longer 
jam dissolution times. 



5.2 Specifying the fundamental diagram directly 

Instead of deriving the fundamental diagram based on the 
approach taken in section |3 we now try to obtain the val- 
ues for the critical densities and capacity flows directly, by 
looking at the STCA's (k,q) diagrams in Fig. (an equiv- 
alent procedure would be to measure the capacity flows 
directly from the time-space diagram in Fig. . Consider- 
ing the STCA's (k,q) diagrams in Fig.r^O we can estimate 
its capacities at approximately qca,p = 0.34 vehicles/time 

step ~ 1220 vehicles/hour, and <^ a p = 0.67 vehicles/time 

step w 2400 vehicles/hour for w max = f and v^ax = 5 
cells / time step, respectively. The stochastic noise p was set 
to 0.1 for all three segments. Changing p to 0.5 for these 
segments, we can estimate the capacities at approximately 
<?cap = 0-15 vehicles/time step « 540 vehicles/hour, and 

9cap = 0.34 vehicles/time step « f220 vehicles/hour for 
Umax = 1 and %ax = 5 cells/time step, respectively. 

Instead of calculating the capacity flows from the av- 
erage free-flow speeds and the critical densities, as was 
done by means of equation (|23|l . we now specify these ca- 
pacity flows directly to the LWR's fundamental diagrams 
and calculate the critical densities from them. The results 
we obtained, are visualized in the time-space diagrams 
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Fig. 11. The (k,q) fundamental diagrams for the STCA 
model. Left: two diagrams for Umax = 1 cells/time step. Right: 
two diagrams for u^^x — 5 cells/time step. Each time, the 
slowdown probability p € {0.1,0.5}. Note how a slower max- 
imum speed makes the diagrams more curved, and how an 
increasing slowdown probability leads to both a lower critical 
density and capacity flow. 

of Fig. [RJ Because the STCA's capacity flows are now 
better approximated, and not overestimated as with the 
previous methodology, there seems to be a better quali- 
tative agreement for both noise levels with the STCA's 
time-space diagrams in the left parts of Fig. and Fig. |H1 




Fig. 12. Time-space diagrams showing the propagation of 
densities during 3000 time steps for the road in the case study. 
Left: the results for the LWR model with p — 0.1. Right: the 
results for the LWR model with p = 0.5. In both diagrams, the 
fundamental diagram was specified directly to the LWR model, 
by explicitly stipulating the capacity flows of the STCA model. 
As a result, there seems to be a better qualitative agreement 
for both noise levels with the STCA's time-space diagrams in 
the left parts of Fig.|7|and Fig. OH 

In Fig. [CH we have depicted the absolute differences 
between this approach and the STCA's time-space dia- 
grams. Comparing this to the previous results of Fig. EH 
we can see that in both cases the buildup and dissolution 
of congestion queues is in good qualitative agreement for 
both noise levels. As such, we come the conclusion that it 
is vital to correctly capture the capacity flows of the STCA 
model. Neglecting this property, can result in severe dis- 
tortion of the system dynamics for higher noise levels. 

Note that the LWR model is able to correctly capture 
the first-order effects of jam buildup and dissolution, and 
that, due to its microscopic treatment, the STCA model 
allows us to visualize the higher-order effects inside jam. 
However, as is evidenced by this and the previous sec- 
tion, it is very important to correctly capture the capacity 
flows in the STCA model, otherwise a growing discrep- 




Fig. 13. Time-space diagrams showing the differences in den- 
sities for the STCA and LWR models, during 3000 time steps 
for the road in the case study. Darker regions indicate large 
differences between both modeling approaches. Left: the dif- 
ferences for p = 0.1. Right: the differences for p = 0.5. In both 
cases, the differences are less pronounced, showing only dark 
edges at the bottom of the triangular-shaped region in segment 
A. 

ancy between the LWR and STCA model is introduced 
with higher noise levels. 



6 Conclusions 

In this paper, we presented an alternate methodology for 
implicitly incorporating the STCA's stochasticity into the 
macroscopic first-order LWR model. The innovative as- 
pect of our approach, is that we derive the LWR's funda- 
mental diagram directly from the STCA's rule set, by as- 
suming a stationarity condition that converts the STCA's 
rules into a set of linear inequalities. In turn, these con- 
straints define the shape of the fundamental diagram that 
is then specified to the LWR model. 

For noise-free systems, our method is exact. In the 
presence of noise, however, the capacity flows in the de- 
rived fundamental diagram are overestimates of those of 
the STCA model. This discrepancy can be explained as 
follows: the underlying assumption for the LWR model, is 
that the fundamental diagram is assumed to be exact, and 
implicitly obeyed, i.e., the existing equilibrium relation is 
representative for the real traffic situation. In the origi- 
nal LWR formulation, this relation was also assumed to 
hold also for non-stationary traffic (which is a more or less 
reasonable assumption if we consider long and crowded 
roads). Our calculations have shown that a direct transla- 
tion of the STCA's rule set into the LWR's fundamental 
diagram, does not always result in a valid fundamental di- 
agram, especially for higher noise levels. As such, there can 
be a significant difference between an average fundamental 
diagram (STCA) and a stationary fundamental diagram 
(LWR). As a result, the STCA model is able to temporar- 
ily operate under larger flows and densities than those 
possible for the LWR's stationary fundamental diagram. 
A logical course of action would to better approximate the 
STCA's fundamental diagram. By doing so however, we 
lose the advantage gained through an explicit derivation 
of the fundamental diagrams' outer envelope, almost cer- 
tainly leading to extra conditions that need to make fur- 
ther assumptions about its shape. Directly specifying the 
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STCA's capacity flows to the triangular LWR fundamen- 
tal diagram, effectively remedies most of the mismatches 
between both STCA and LWR models. 

Our methodology sees the STCA complementary to 
the LWR model and vice versa, so the results can be 
of great assistance when interpreting the traffic dynam- 
ics in both models. Especially appealing, is the fact that 
the STCA can visualise the higher-order characteristics of 
traffic stream dynamics, e.g., the fans of rarefaction waves. 
Nevertheless, because the LWR model is only a coarse rep- 
resentation of reality, there are still some mismatches be- 
tween the two approaches. One of the main concerns the 
authors discovered, is as hinted at earlier, the fact that us- 
ing a stationary fundamental diagram (i.e., an equilibrium 
relation between density and flow), always overestimates 
the practical capacity of a stochastic cellular automaton 
model. As such, it is vital to correctly capture the capac- 
ity flows in both STCA and LWR models, a remark that 
we feel is valid for all case studies. 
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